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An analysis of the Dalitz plots for the decays D"*" K~ti'^-k^ and n~n^n^ indicates 

that structure with significant phase variation is required in the s-wave K~-k^ and tt'tt^ 
systems in the region below 1 GeV/c^ effective mass. A constant, non resonant "contact" 
term, together with a simple s-wave Breit Wigner amplitude with M„ = 797 ± f 9 ± 43, r„ = 
175±12±f2 MeV/c^ for the K--K+ system and = 478123 ±17, r„ = 324tll±21 MeV/c^ 
for the TT^Tv^ system provide good fits to the data in this region. 



1 Introduction 

Decays of D mesons to three pseudo scalar mesons can produce intermediate, two meson systems 
with natural parity. Kinematics and angular momentum barrier factors generally favour scalar 
(J^ = O"*") over vector 1~ or tensor 2^ systems. Analysis of these decays can therefore provide 
new information on the poorly understood scalar meson system. 

In most instances, D decays appear to proceed through quasi two body channels. In the cases 
of decays ^ to 7r"_7r+7r~'" and to K~7r~^7r~^ however, intermediate states found in the Particle 
Data Group listings a (PDG) have been found to contribute only partially to the decays, with 
three body modes apparently dominant!' 1 We report here on analysesill of these decays using 
the largest samples so far examined. 
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2 Fermilab Experiment E791 

In this experimenti 2x lO^*' charm enhanced interactions between 500 GeV/c vr^ mesons and thin 
foil targets (four C, one Pt) are recordedQ. Foils are spaced so that decays of D mesons produced 
in them occur predominantly between them. Silicon strip detectors measure D decay lengths 
(mean ~5 mm) with a typical precision of 350 fim. Charged particle momenta are measured 
with two dipole magnets and drift chambers and two threshold gas Cherenkov detectors provide 
K/ir/p separation in the 6-60 GeV/c momentum range. 

Events with one negative and two positive charged tracks making a good vertex, well sep- 
arated from any foil and from the primary interaction point, are selected for further analysis. 
The resultant momentum is required to point back to within 40 /im of the primary vertex. For 
the i^~7r+7r+ sample the K~ is required to have momentum 6 GeV/c< \pk\ < 40 GeV/c and 
to be identified by the Cherenkov detectors. Other cuts are also made to reduce contamination 
from tracks that are really associated with the primary vertex. With these selection criteria, 
the major sources of background are mostly misidentified (or incomplete) charm decays, with 
a small, "combinatorial" contribution from false secondary vertices. Yields of 15,090 K'^tt^tt^ 
(6% background) and 1,686 tt^tt^tt^ (28% background) remain in the signal regions. 



3 Isobar Analysis of D Decay to Three Pseudo-Scalars 

Decays of D mesons to three pseudo scalars i, j, k can be described approximately as a sum 
of isobar amplitudes, each corresponding to a quasi two body decay D — > i?(— > i j) k. Each 
amplitude must satisfy Lorentz invariance and conserve total spin and has the form 



F^{q,r^)F^{p,rj,) x BWn.{s^,) x {-2y \p\^\q]-^ Pj {cos 



where p, q are momenta of i and k respectively in the {ij) rest frame and cos = p ■ q. Form 
factors -Fo (-Fr) for D (R) are parametrized in terms of effective radii Ro (Rr)- A Breit Wigner 
propagator BW^ = [s^ — Sij — i^/s]^ T{sij)] ^ describes the resonance R with spin J, mass ^/s^. 
Suffix R denotes a quantity evaluated at Sij = s^. 

A decay is kinematically specified by coordinates {sij,Sik) on a Dalitz plot and by the 
three body mass Mj^-fc whose distribution about the D mass mu due to resolution is G{Mijk), 
determined from the data. The signal distribution of decays is 



'Ps{sij,Sik,Mijk) = G{Mijk) X 



ikj 



in which three body decays (NR) are described by a constant, "contact" amplitude Ojvjje*'^'^^ . 

The distribution VBisij, Sik, Mijk) of backgrounds is modelled using events from side bands 
in Mjjfc for combinatorial, and from MC simulations for broken charm sources. The efficiency 
e{sij,Sik) for observing decays is determined from MC simulations. 

Complex coefficients ae^^ and resonance parameters are determined from fits to a linear 
combination of ePs and Vb in the correct proportions. The fraction of D decay to an isobar R 



is fn = G{Mijk) J {Anl^dsijdsik/ J Vgdsijdsik evaluated at Mijk 



m, 



4 The K-TT+-K+ Channel. 



The Dalitz plot in figure ^(a) shows clear bands in the iiri(890) regions with an asymmetric 
structure indicative of interference with a significant, underlying s-wave. An accumulation of 
events near 2 (GeV/c^)^ indicates contributions from i^^(1430) or /C2(1430). A fit is made 
including isobar amplitudes for these resonances and the i^i(1688), each with mass and width 
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Figure 1: (a) Dalitz plot for — > K'tt^tt'^ decays. Projections for (b) smaller and (c) larger Kn effective mass 
for data (crosses) and fit (line) to model I described in the text. Arrows indicate regions where the fit is worst. 

fixed at their PDG values. Results are shown in table |^ (model I) and compared with the K~Tr~^ 
mass projections in figures Hb) and (c). As in previous analyses afl using this model, a large 
NR component is required (~ 90%) and the sum of all contributions is ~ 140% indicating much 
destructive interference. The fit is poor with a x^i found by comparing predictions of model I 



Table 1: Resonant fractions / for D K tt tt , and tt tt tt . Models I-IV are described in the text. 



D+ ^ 


K 7r+7r+ 




ifi(890)7r 


i^o(1430)7r 


if2(1430)7r 


A'i(1688)7r 


I 


90.0 ±2.6 




13.8 ±0.5 


30.6 ± 1.6 


0.4 ±0.1 


3.2 ±0.3 


II 


13.0 ±5.8 


47.8 ± 12.1 


12.3 ± 1.0 


12.5 ± 1.4 


0.5 ±0.1 


2.5 ±0.7 




±4.4 


±5.3 


±0.9 


±0.5 


±0.2 


±0.3 


D+ ^ 


7r~7r"^7r^ 


CTTT 


p(770)7r 


/o(980)^ 


/2(1270)^ 


p(1450)7r 


III 


38.6 ± 1.4 




20.8 ±2.3 


7.4 ±4.3 


6.3 ±3.3 


22.6 ±2.1 


IV 


7.8 ± 6.0 


46.3 ± 9.0 


33.6 ±3.2 


6.2 ± 1.3 


19.4 ± 2.5 


0.7 ±0.7 




±2.7 


±2.1 


±2.2 


±0.4 


±0.4 


±0.3 



with data in equal square bins in the Dalitz plot, of 167 for 63 degrees of freedom u. Major 
discrepancies are at low K~tt^ masses, arid the corresponding high mass reflection, and can be 
seen in earlier data for both this channels and for D° — > K~tt^'k°u. However the large E791 
sample underlines their significance, and demands that a better model be found. 

First, the radii Rr,, Rji and the mass and width of -ftro(1430) are allowed to float This 
produces little improvement. Next, a further scalar amplitude "k" having a Breit Wigner form 
similar to Kq{1430) is added (model II in Table ||). This does provide an excellent description 
of the data with = 46/63. The k mode is dominant while the NR component is consistent 
with zero. The sum of fractions is ~ 86%. The k mass and width, determined from the fit, are 
M„ = 797 ± 19 ± 43 and r„ = 175 ± 12 ± 12 MeV/c^. In this model, however, the Ko{U30) 
mass and width are found to be 1459 ± 7 ± 12 MeV/c^ and 175 ± 12 ± 12 MeV/c^ respectively, 
differing somewhat from the values from LASsi. 

The s-wave BW form used for the k in model II appears to be important. Variations in 
this, or in form factors for k and NR terms, are found to affect principally the amounts of k 
and NR, but to have relatively small effects on k mass and width parameters. Replacing the 
complex BW in model II by its absolute magnitude (a "k" mass enhancement, but no phase 
variation) results in unphysically large k and NR fractions. A vector form for k converges with 
mass ~ 1100 MeV/c^ and negligible fraction. A tensor k fit fails to converge. 



5 —>■ vr 7r+7r+ Decays. 

Figure |2|(a) is the Dahtz plot for these decays @. Table || shows results for model III in which 
only resonances from the PDG are included. The fit is poor (x^/^ = 1-5), especially in the low 
mass TT~7r~^ system as seen in figure |2|(b). The contribution from p(1450) is larger than that 
from p(770) - a somewhat unlikely result. As observed in the previous analysistJ of this channel 
with a similar model, the NR mode dominates. 
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Figure 2: (a) Dalitz plot for D'^ 7r~7r+7r"'" decays, (b) Projection of ■K~Tr'^ effective mass. Data (points with 
errors) and predictions (solid line) from model III described in the text indicate poor agreement at low mass. 

Adding a scalar amplitude {a, similar to the k in model II) results in a superior fit jv = 
0.9). Results for this (model IV) are in table |l[ Mass and width for a are found to be = 
478t23±17 and = 324l4o±21 MeV/c^. Both NR and p(1450)7r contributions are insignificant, 
and the an amplitude is dominant. As for the k, the s-wave BW form seems to be the best 
model of those applied to these data. 



6 Conclusions. 



We conclude that the s-wave K~tt~^ system from decay of to iC~7r^7r+ has amplitude and 
phase variation consistent with a sum of NR, ii'o(1430) and scalar k terms. Similarly, the s-wave 
7r^7r+ system from decay of vr^Tr+vr"'" has amplitude and phase variation consistent with 

a sum of NR, /o(980) and scalar a terms. Though this hints at the existence of scalar k and a 
resonances, use of Breit Wigner functions to describe such broad states near threshold may not 
be appropriate, and other interpretations of the data may be possible. 
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